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The interaction between doped graphene nanoislands connected by narrow junctions constitutes 
an ideal testbed to probe quantum effects in plasmonic systems. Here, the interaction between 
graphene plasmons in neighboring nanoislands is predicted to be extremely sensitive to the size 
and shape of the junctions. The reported ab initio calculations reveal three different regimes of 
interaction: (1) for narrow bridges (< 4 carbon-atom rows), the conductance of the junction is too 
low to allow electron transport and the optical response is dominated by a characteristic bonding 
dipolar dimer mode that also appears in a classical description; (2) for wider junctions (4-8 carbon 
rows), a strong charge polarization is induced across the junction, which gives rise to a novel junction 
plasmon that has no counterpart in a classical description; (3) for even wider junctions (> 8 rows), 
their conductance is sufficiently large to allow charge transport between the two graphene islands, 
resulting in a pronounced charge-transfer plasmon, which can also be described classically. This 
work opens a new path for the investigation of intrinsic plasmon quantum effects. 



I. INTRODUCTION 

A major advantage of plasmons -the collective exci- 
tations of conduction electrons in metals- is that they 
can produce high levels of light confìnement and optical 
field enhancement |T], which enable them to strongly in- 
teract with objects that have relatively low optical cross- 
sections, such as atoms and simple molecules in solid- 
state environments 0. Additionally, confined plasmons 
can display intense dipoles through which they couple to 
external radiation. Therefore, plasmons can act as inter- 
mediaries (antennas) that boost the interaction of light 
with molecules (e.g., in techniques as useful as surface- 
enhance Raman scattering (SERS) [31 [4] and surface- 
enhanced infrared resonant absorption (SEIRA) [5], and 
also in promising applications to tumor removal [6l , drug 
delivery jjj, improved photovoltaics [3, and catalysis [9 ). 
An interesting scenario is found when a plasmon is sup- 
ported by a single molecular structure rather than an ex- 
tended nanoparticle. Many molecules have been found to 
exhibit plasmons p!QHl4] , including the Geo molecule and 
carbon nanotubes [151 [16]. Like other carbon allotropes, 
they display plasmon bands at energies around 5eV 
(tt plasmons) and ^ 15 — 20 eV (a plasmons). However, 
these are conventional plasmons, in the sense that they 
are not too sensitive to the charge state of the molecules, 
similar to those found in noble metal nanoparticles. In- 
stead, we are here concerned with plasmons hosted by 
doped graphene molecular nanostructures, which in con- 
trast to noble-metal plasmons, only occur when they are 
electrically charged. These plasmons emerge at lower en- 
ergies for typical levels of doping. 

Graphene has recently emerged as an attractive elec- 
trically tunable, optical material [17 displaying strong 
infrared plasmons, the frequencies of which scale roughly 
as |n|^/^ with the doping charge density n pISHST] . 
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Graphene plasmons provide unprecedented levels of light 
confìnement and fìeld enhancement [22], which are 
promising tools for accessing quantum-optics phenom- 
ena [23H26] . In this context, graphene offers the addi- 
tional benefìts of a robust structure and a simple and 
convenient electrical tunability for controlling quantum 
interactions [27 . The electrical modulation of graphene- 
plasmon-related phenomena also suggests potential ap- 
plication in optical waveguiding and signal processing 
[28-30 , metamaterials [31^ 32 , quantum optics [22 , and 
spectrometry [33] . Experimental evidence of electrostatic 
control over plasmonic features in the absorption spectra 
of graphene has been recently reported [33j, followed by 
the observation and spatial mapping of confined plas- 
mons in this material [34, 35j. Additionaly, graphene 
has been used to extrinsically modulate the plasmons of 
neighboring metallic structures [36^, '37| . These advances 
have sparked considerable interest into the unique optical 
behavior of nanostructured graphene. 

Here, we describe intrinsic quantum effects in the plas- 
monic response of bowtie graphene nanodimers bridged 
by thin junctions. Due to the two- dimensionai character 
of this material, the addition of a small number of atoms 
(< 10) is suffìcient to dramatically modify the absorption 
spectrum of the entire dimer. Three distinct regimes of 
plasmonic behavior are predicted by our fìrst-principles 
calculations. For separated nanotriangles bridged by nar- 
row junctions (< 4 atomic rows wide), we find a charac- 
teristic hybridized bonding dipolar plasmon (BDP) mode 
[38] , In the limit of wide junctions (more than 8 atomic 
rows), a dramatic reshaping of the optical spectrum leads 
to a pronounced charge-transfer plasmon (GTP), which 
emerges at lower energies, similar to what happens when 
the gap vanishes in a conventional metallic nanoparti- 
cle dimer [39-42^. The plasmonic features in both of 
these regimes are also found within a classical electrody- 
namics description using the appropriate dielectric per- 
mittivity for the graphene, although the behavior of the 
GTP is substantially modifìed by quantum mechanical ef- 
fects (see below). In contrast, for intermediate junction 
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FIG. 1: Quantum junction plasmons in graphene bowties. (a) Details of the junction region in the structures under 
study, with definitions of bridge length n and width m. (b-d) Spectra for selected bridge widths m = 0, 4, 8, and various 
lengths as indicated by different colors. The insets show the complete graphene structures for n = 2. (e) Plasmon resonances 
as a function of bridge width m. The color code for different lengths n is given in the upper inset. The area of the circles is 
proportional to the area under the extinction peak for each plasmon feature. Classical-electrodynamics plasmon energies are 
shown by dashed curves for n = 4 (red) and n = 8 (blue), flanked by shaded areas that represent the strength of the modes. 
(f) Plasmon resonances as a function of bridge length n. The graphene is doped to a Fermi energy Ef = 0.4 eV and has a 
mobility /x = 10, 000cm^/(V s). The full length of the bowties is 8nm. 



widths (4 — 8 atomic rows), a plasmonic feature shows 
up at intermediate energies, which is absent in classical- 
electrodynamics calculations. We refer to this mode as a 
junction plasmon (JP) because it is caused by quantum 
effects in the small molecular junction. We also predict 
a minor dependence on the length of the junction and 
a strong dependence on its width, thus confìrming the 
importance of a quantum mechanical description of the 
optical properties of bridged graphene dimers. 



II. RESULTS AND DISCUSSION 

We simulate the optical response of nanoscale bowtie 
graphene structures using the random-phase approxima- 
tion (RPA) and wave functions derived from a tight- 
binding description of the carbon sheet. Further details 
of the method are given elsewhere [43]. The triangular 
structures forming the bowties are oriented with respect 
to the graphene lattice in such a way that they have 
armchair edges, which prevent undesired losses typically 
observed in zigzag-edge structures, and caused by the 
presence of zero-energy electronic edge states [43 . Such 
States are not present in the structures under considera- 
tion, because they have the same number of atoms in 
both carbon sublattices [44 . For comparison, we in- 
clude classical-electrodynamics calculations, which are 
performed by describing the graphene through its local- 
RPA surface conductivity ^22^ using a finite elements 



method (COMSOL). We fìx the side length of the nano- 
triangles to 8 nm in ali cases (~ 10^ atoms in each trian- 
gle). We set the Fermi energy of our structures to Ep = 
0.4 eV, and the intrinsic damping to hr~^ = 1.6 meV, 
corresponding to a DC mobility of 10,000 cm^/(Vs) [21]. 

Our main results are summarized in Fig. [ij We show 
details of the junction region in Fig. [TJa) for three char- 
acteristic values of the width (m = 0, 4, 8, where m is the 
number of carbon- atom zigzag rows forming the junction, 
and m = stands for the non-touching confìguration) . 
The corresponding spectra are shown in Fig. [TJb)-(d) for 
a range of junction lengths n = 2 — 10 [n is defìned as the 
number of carbon hexagons that are needed to join the 
graphene triangles for m = 2, see Fig.[l];a)]. These spec- 
tra show a clear trend from high to low energy features as 
the junction width is increased, irrespective of the length. 
This behavior is further examined in Fig.jlje), where the 
plasmon features are arranged as a function of energy 
(vertical axis) and junction width (m, horizontal axis) 
for different values of the junction length n (see color 
code in the upper inset; full spectra are provided in the 
Appendix (Fig. |4|. The area of the circles is made pro- 
portional to the intensity of the plasmon mode defìned as 
the area under the corresponding plasmon peaks in the 
extinction spectra (for light polarization along the gap). 
For narrow bridges, the spectra exhibit prominent plas- 
monic features centered around ~ 0.47 eV. In contrast, 
the spectra for the wide-junction limit are dominated 
by lower-energy plasmonic features around ~ 0.22 eV. 
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FIG. 2: Plasmon induced- charge distributions for representative graphene bowties. (a-d) Density plots with the 
color of each atom indicating its induced charge for difFerent junction widths m = 0,2,6,8 and the same length n = 4. (e) 
Induced charge integrated along the horizontal direction and averaged over four nearest carbon-atom neighbors. 



At intermediate values of the junction width, there is 
a complex transition between these two regimes, with 
intermediate-energy features around ~ 0.35 eV showing 
up in the spectra. The transition is fast, but not singu- 
lar. These conclusions seem to be rather independent of 
the length of the junctions n. As Fig. [TJe) illustrates, no 
signifìcant spectral changes are induced when the plas- 
mon energies are plotted as a function of junction length 
n. The three distinct behaviors of the plasmon energies 
noted above for narrow, intermediate, and wide junctions 
appear for ali values of the junction length. 

Further insight into the character of these plasmons 
is provided by examining their induced-charge distribu- 
tions, which are plotted in Fig. [2|a)-(d) for representa- 
tive structures along the transition from narrow to wide 
junction regimes. The polarization profìle for the high- 
energy plasmon around ~ 0.47 eV in the non-touching 
dimer [m = 0, Fig. [2|a)] displays a dipole-dipole pat- 
tern that is not substantially modifìed when a narrow 
junction [m = 2, Fig. [2|b)] is present between the two 
triangular islands of the structure. This distinct polar- 
ization pattern is characteristic of a BDP. In contrast, 
the intermediate energy plasmon around ~ 0.35 eV [Fig. 
[TJc)] for a junction width m = 6 shows a clear dipolar 
polarization pattern across the junction. This is the JP 
that is formed from the locai electronic properties of the 
junction. A weak charge transfer between the two nan- 
otriangles is also evident in the plot and corresponds to 
the CTP, which is delocalized over the entire graphene 
structure like an standing wave. The latter shows a con- 
sistent H — charge polarization that does not change with 
increasing m, as shown for m = 8 in Fig. [2|d). The CTP 
becomes the dominant feature once the junction is suf- 
ficiently wide. The signs superimposed on the density 
plots of Fig. |2|a)-(d) qualitatively correspond to the dis- 



tribution of the plasmon induced charge as a function 
distance to the dimer center, which is shown in Fig. [2|e) 
(see also Fig. [6] in the Appendix). This analysis is rather 
independent of the length of the junction n. An in-depth 
inspection of bowties of increasing length leads to results 
that are consistent with the above conclusions, including 
the emergence of quantum junction plasmons once the 
width of the junction and the length of the structure are 
suffìciently large (see Fig. 10 in the Appendix for more 
details). 

A classical-electrodynamics description of these 
graphene nanostructures (see broken curves in Fig. [TJe) 
and Fig. [s] in the Appendix) does not completely repro- 
duce the plasmonic behavior derived from fìrst-principles. 
For instance, classical calculations miss the intermediate- 
energy junction plasmon that we observe in the quantum 
calculations for intermediate junction widths. Another 
discrepancy is that classical theory predicts a smooth 
fadeout of the dipole-dipole mode exhibited by non- 
overlapping triangles when the junction width is in- 
creased, while the CTP shows a singular behavior, as 
it migrates towards zero energy in the limit of vanishing 
junction width. This is in striking contrast to the re- 
sults from the quantum description, in which the CTP is 
rather dispersionless and fades out smoothly in that limit 
[see Fig. [TJe)]. The dipole-dipole and charge-transfer na- 
ture of the classically predicted plasmons is clear from 
their induced-charge distributions (see Appendix, Fig.[5|. 

The electronic states involved in the plasmon of sep- 
arated nanotriangles [Fig. [sj^a)] are nearly identical 
to those of individuai triangular islands, but with a 
small amount of Coulombic interaction and hybridiza- 
tion. When a bridge of intermediate width is added to 
form a junction [Fig. jsjb)], the level of hybridization of 
these states increases, and new electronic junction states 
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FIG. 3: Electronic junction states and electron- hole transition strengths. (a,b) Electron density distribution of 
electronic states in n = 4 bowtie antennas for m = ((a), separated triangles) and m = 6 ((b), bridged triangles). For each 
state of energy E > 0, there is a state of energy —E (not shown) with identical density distribution. Only low-order states 
dose to the Dirac point (zero energy) are shown. The energies of these states are indicated by black lines right under the 
corresponding density plots. (c) Dipole matrix elements between electronic states of the same bowtie as in (b). The area of the 
circles is proportional to the dipole strength. We only show transitions between states separated by energy differences larger 
than \Ei — Ej\ > 0.01 eV. The solid lines connect initial and final energies corresponding to 0.47 eV (BDP), 0.35 eV (JP), and 
0.22 eV (CTP) excitations. 



emerge. In particular, two new junction states are ob- 
served near the Dirac point (zero energy), which are ex- 
pected due to the presence of carbon zigzag edges in the 
bridge (45] [46]. Thus, the junction plasmon noted above 
must be supported by excitations involving electrons or 
holes in these electronic junction states. The strength of 
the electron-hole pair (e-h) dipole transitions in this sys- 
tem are represented in Fig. [sjc) as a function of initial 
and final energies (the strength is represented through 
the area of the circles). Clearly, the plasmon energies 
(solid curves) do not strongly overlap with the dominant 
e-h transitions. This energy mismatch shows that the op- 
tical transitions are not single-particle excitations and is 
a strong argument supporting our conclusion of collective 
plasmonic nature of the bowtie optical excitations under 
consideration. Similar results are observed for different 
widths of the junction (see Appendix, Fig. [9|. 

III. OUTLOOK AND CONCLUSIONS 

Graphene dimers provide advantageous systems com- 
pared to metallic dimers both because of the struc- 
tural robustness of this material and because the car- 
bon structure can be reliably imaged with atomic de- 
tail [47l|48], and therefore, a detailed correlation between 
atomic structure and optical response is within reach ex- 
perimentally. This is not the case of gaps formed by 
closely spaced three-dimentional metallic particles or tips 
[49^ i50j , which suffer from inherent faceting and struc- 



tural uncertainties. Planar graphene nanostructures are 
thus ideal systems on which to study the interplay be- 
tween structural features (e.g., ripples and defects) and 
the plasmonic response (for example, the amount by 
which they increase the plasmon scattering and dephas- 
ing rate via electron-hole pair generation). Recent ad- 
vances in atomically resolved imaging [ST and patterning 
[52 of graphene should facilitate this task. 

The remarkably large variations here reported in the 
optical properties of plasmonic structures consisting of 
thousands of carbon atoms by just adding a few (< 10) 
extra atoms provides a handle for tuning plasmons in 
graphene. Together with the already established chemi- 
cal [53] and electrical [54 tunability of this material, this 
adds a novel tuning parameter, which can be potentially 
exploited for ultrasensitive optical sensing (e.g., by ana- 
lyzing the dependence on the presence of molecules near 
the junction, which can modify the locai valence charge 
density). 

Although for simplicity we only discuss above free- 
standing graphene with a Fermi energy Ep = 0.4 eV, 
no qualitative changes are observed for different doping 
or when the graphene sheet is placed on a substrate (see 
Appendix, Figs. [7|and[8|. Quantum effects are found to 
be more significant when the level of doping is reduced 
(see results for Ep = 0.2 eV in the Appendix, Figs.[7|), as 
one expects for the smaller number of doping electrons 
participating in the plasmons, whereas the CTP energies 
are closer to the classical description when the doping 
increases (see Ep = 0.8 eV results in the Appendix, Figs. 
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In conclusion, using a fully quantum mechanical ap- 
proach we have studied the optical properties of graphene 
dimers consisting of several thousand carbon atoms con- 
nected by a narrow graphene junction. The optical prop- 
erties of the system are found to be strongly sensitive 
to the detailed structure of the junction, with the spec- 
trum exhibiting drastic changes with the addition of a 
few carbon atoms in that region. In addition to the con- 
ventional hybridized and charge transfer plasmons ob- 
served for metahic dimers [39-42 , graphene dimers dis- 
play junction plasmons of a purely quantum mechanical 
origin, which does not seem to appear in conventional 
metallic junctions [55]. The plasmonic response of this 
graphene structure is found to be strongly influenced by 
quantum effects, making it of signifìcant interest in the 
emerging fìeld of quantum plasmonics as well as for novel 
nanophotonics and optoelectronics applications. 
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APPENDIX 

In this appendix, we provide detailed results of 
first-principles extinction spectra, electron states, and 
plasmon-induced charge densities for the graphene 
bowtie structures under consideration. Additionally, we 
study the evolution of plasmons in bowties of increasing 
length. Classical-electrodynamics calculations are also 
offered for comparison. 

We choose a Fermi energy Ep = 0.4 eV (except in Fig. 
[t]) and an intrinsic damping Tìt~^ = 1.6 meV, correspond- 
ing to a DC mobility of 10,000 cm^/ (Vs). We set the tem- 
perature to T = 300 K, although thermal effects are neg- 
ligible for the plasmons under consideration, the energies 
of which are high compared with ksT. The external light 
polarization is taken along the junction between the nan- 
otriangles forming the bowties. The triangle side- length 
is set to 8nm, except in Fig. [lO] Our fìrst-principles 
calculations are performed following methods described 
elsewhere [43]. For the classical calculations of Fig. [sjwe 
use a fìnite-element method (COMSOL) following a pre- 
viously reported prescription [22 according to which the 
graphene is treated as a thin dielectric layer. 

We show in Fig. [4] the full series of extinction spectra 
from which Figs.[TJe) andjljf) are constructed. The cross 
section is given in square nanometers, so it can be direct ly 
compared with the geometrical area of the graphene in 



the bowties ~ lllnm^. The maximum cross-section 
obtained in these first-principles calculations is approxi- 
mately half the geometrical area. In contrast, the cross 
sections derived from classical electrodynamics (Fig. [5| 
can reach substantially large values compared with the 
geometrical area. The energies of the plasmons derived 
from classical calculations are represented in Fig. [TJe) for 
comparison. 

The charge densities associated with the plasmons 
(Fig. [2| are examined in more detail in Fig. [6| where 
they are shown to exhibit strong modulations between 
nearest-neighbor carbon-atom rows. These modulations 
oscillate between three smoothly varying limits, as ex- 
pected from the three different types of atomic sites in 
the armchair edges of the nanotriangles under consider- 
ation. The curves plotted in Fig. [2j and reproduced in 
Fig. |6] for convenience, correspond to averages over four 
carbon-atom neighbors. 

We have considered the effect of doping in Fig.[7| where 
we plot data similar to Fig. [TJe) {Ep = 0.4 eV), but 
changing the doping level to Ep = 0.2 eV (Fig. [7|^a)) 
and Ef = 0.8 eV (Fig. [7|^b)). As a general trend, we 
observe a reduction in the role of quantum effects as the 
doping level is increased. For example, for Ef = 0.2 eV 
we observe that charge-transfer plasmons (CTPs) pile 
up around 0.15 eV in the quantum description, whereas 
classical theory places them at much lower energies; in 
contrast, both classical and quantum descriptions agree 
much better in the energies of CTPs for Ef = 0.8 eV; 
and an intermediate situation is encountered for Ef = 
0.4 eV. Likewise, junction plasmons (JPs) appear to lose 
strength relative to the rest of the features as in- 
creases. Actually, these plasmons cover a wide spec- 
tral region extending above the energies of CTPs for 
Ef = 0.2 eV as the junction width is varied, whereas 
they are only discernable in the m = 3 — 7 range of junc- 
tion widths for Ef = 0.8 eV. 

For simplicity, we present in this paper results for self- 
standing graphene, except in Fig. [8| where we study the 
effect of depositing the carbon layer on a glass substrate 
(e = 2). In a classical description of graphene, the net 
effect of the substrate is to change the conductivity a to 
an effective value 2a/{e-\-l) [30J. This leads to a redshift 
of the plasmon energies by a factor - ^^2/(6 + 1). We 
show in Fig. [Sj that a similar trend is observed in small 
graphene structures within our quantum-mechanical de- 
scription, thus corroborating that the conclusions of our 
study also apply to the more physical situation of sup- 
ported graphene nanostructures. The features of the ex- 
tinction spectra are maintained by introducing a sub- 
strate, which essentially produces a redshift that is sub- 
stantially smaller than that predicted by classical the- 
ory in our relatively small structures, particularly for the 
lowest-energy features and the widest junction considered 
in Fig. [8| 

The electron wave functions and the strength of dipo- 
lar transitions connecting them are considered in Fig. [9| 
which is an extension of Fig. [3| including more junction 
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FIG. 4: Extinction spectra in bowties of increasing junction width. We bere show the full extinction spectra from 
which Figs.[l{e) and[l];f) are constructed as a function of junction width m and length n. 



widths. One can clearly observe the increasing involve- 
ment of electron junction states as the junction becomes 
wider. 



Finally, we study in Fig. 10 the evolution of plasmons 
in bowties as a function of the length of the structures. 
We fix the junction length to n = 4 hexagons and con- 
sider three different junction widths consisting of m = 2, 
6, and 8 carbon zigzag rows. In general, we observe a 
migration of modes with similar symmetry towards lower 
plasmon energies as the bowtie length increases. For the 
two shorter series under consideration (half length = 2, 
3, and 4nm), we fìnd a dominant CTP, with opposite 
induced charges in each half of the dimer. For larger 



bowties, higher-order modes become relevant. For exam- 
ple, the structure with the narrower junction (m = 2) ex- 
hibits a dipole-dipole plasmon, which ends up at ~ 0.4 eV 
at the full length 2X ^ 14 nm (lower row of Fig.[lO|. For 
this junction width, the CTP becomes relatively weak 
with X > 6 nm, so it does not make a relevant contribu- 
tion to the full structure. In contrast, the CTP remains 
a dominant feature for ali the lengths under considera- 
tion in structures bridged by wider junctions (m = 6 and 
8). As the length increases, a new dipole mode starts 
developing, which is associated with large charge pileup 
at the junction. This is the intermediate-energy feature 
observed in Fig. [TJe). 
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FIG. 5: Classical-electrodynamic calculations. (a,b) Extinction spectra of bowties for two different values of the junction 
length n (we use the same notation as in Sec.[n]for m and n). (c) Induced charge distributions associated with representative 
plasmons in these structures. Upper row: charge-transfer plasmons. Lower row: dipole-dipole hybridized plasmons. 
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FIG. 6: Plasmon charge density along the bowties. Induced charge associated with the plasmons considered in Fig. [2] 
Dots: Charge density integrated across the width of the bowtie as a function of the position of the carbon atoms relative to 
the center of the junction. Curves: Same, averaged over four nearest-neighbor atom positions along the axis of the bowtie. 
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FIG. 7: EfFect of doping. We show the plasmon energies and strengths for two difFerent levels of doping (Ef = 0.2 eV in (a) 
and Ef — 0.8 eV in (b)) with the same notation and style as in Fig. [T]^e). The insets show plasmon charge distributions for 
selected geometries. 
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FIG. 9: Electron states and electron-hole transition strengths. We show bere a more extended version of Fig. [3] with 
electron states and electron-hole transition strengths plotted for various values of the junction width m and length n. We only 
show transitions between states separated by energy differences larger than \Ei — Ej\ > 0.01 eV. 
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FIG. 10: Formation of difFerent types of plasmons in bowtie antennas of increasing length. Left: extinction spectra 
for three different values of the junction width m and increasingly larger bowties. The half-length of the bowtie is denoted X 
(see upper inset). The full structure corresponds to a triangle side- length of 8nm (i.e., X ^ 7nm). Right: plasmon charge 
density distributions for the features indicated by arrows in the extinction spectra (in the same order from left to right within 
each row). 



